Unusual quasi-one-dimensional electron dispersions in the spin-1/2 quantum magnet 
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We have studied the electronic structure of the spin-1/2 quantum magnet TiOCl by polarization- 
dependent momentum-resolved photoelectron spectroscopy. From that, we confirm the quasi-one- 
dimensional nature of the electronic structure along the crystallographic b-axis and find no evidence 
for sizable phonon-induced orbital fluctuations as origin for the non-canonical phenomenology of 
the spin-Peierls transition in this compound. A comparison of the experimental data to our own 
LDA+U and Hubbard model calculations reveals a striking lack of understanding regarding the 
quasi-one-dimensional electron dispersions in the normal state of this compound. 

PACS numbers: 71.30.+h,71.27.+a,79.60.-i 



I. INTRODUCTION 



Strongly correlated low-dimensional quantum spin-1/2 
systems have attracted much attention during the last 
years since they can serve as a testing ground for theo- 
retical concepts which are believed to play a key role 
in understanding materials of paramount technological 
potential like the high-T c superconductors or the CMR 
(colossal magnetoresistance) manganitesiSi This applies 
in particular to one-dimensional systems for which ex- 
actly solvable many-body models exists Meanwhile the 
narrowed view on the isolated subsystems, e.g. the elec- 
tron system, has given way to a synopsis of the complex 
interplay between all relevant degrees of freedom, i.e. 
charge, spin, orbitals, and lattice. In this perspective, 
the layered Mott insulator TiOCl has recently aroused a 
lot of interest. 

TiOCl crystallizes in an orthorhombic quasi-two- 
dimensional structure of the FeOCl type, where buckled 
bilayers of Ti-0 are separated by CI ions^ The bilayers, 
which are stacked along the crystallographic c-axis, only 
weakly interact through van der Waals forces. In a local 
picture, the electronic structure is determined by the oc- 
tahedral coordination of the Ti ion in a 3d 1 configuration. 
The strongly distorted octahedra are formed by four O 
and two CI ions, and share corners along the a-axis and 
edges along the 6-axis. Thus, the low-lying charge excita- 
tions occur within the Ti 3d t2 g triplet. This simple pic- 
ture is essentially confirmed by LDA+U calculations, 5,6 
which identify the d xy derived band as slightly split off 
from the bands with d xz and d yz character. 7 Both the 
LDA+U results as well as the observation of a Bonner- 
Fishcr-type magnetic susceptibility at high temperatures 
indicate the existence of spin-1/2 Heisenberg chains in 
this material, mediated by direct exchange along the b- 



axis (J k 660 K), 5 Correspondingly, the sudden drop of 
the susceptibility at T Cl = 67 K to almost zero and a kink 
anomaly at T C2 = 91 K have been discussed in terms of 
an unusual first-order spin-Peierls transition and possibly 
some kind of precursor transition, respectively. While 
the dimerized nature of the ground state indeed could be 
proven by x-ray diffraction^ the results from magnetic 
resonance Ramanii and infrared spectroscopy^ and 
specific heat measurements^ point to the importance 
of strong spin and/or orbital fluctuations in the high- 
temperature phase up to 130 K. In partial contradiction, 
there is recent evidence from cluster calculations in con- 
nection with polarization-dependent optical data that the 
orbital degrees of freedom are actually quenched^ Based 
on Ginzburg-Landau arguments it was concluded that 
the frustration of the interchain interactions in the bi- 
layers give rise to incommensurate order in the interme- 
diate phase, which commensurately locks in below T Cl . 
Such behavior has indeed been observed in the related 
compound TiOBr. 1 - Still, a complete understanding of 
the relevant competing interactions in this compound 
demands for a thorough investigation of the electronic 
properties in the normal state. This seems particularly 
interesting under the perspective of doping charge car- 
riers into the system to drive the chains metallic and 
possibly even superconducting. 16 ' 17 

In this paper we focus on the normal state elec- 
tronic properties of TiOCl at room-temperature and 
slightly above, i.e. well above the transition tempera- 
tures T Cl and T C2 . To this end we performed polariza- 
tion dependent angle-resolved photoemission measure- 
ments (ARPES) which provide momentum-resolved in- 
formation on the electronic structure. The experiments 
are complemented by LDA+U and Hubbard model cal- 
culations. Our photoemission data confirm the one- 
dimensionality of the electronic structure and give no evi- 
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FIG. 1: (Color online) Magnetic susceptibility of TiOCl mea- 
sured by SQUID magnetometry. 



dence for any significant (phonon-induccd) admixtures of 
d xz and d yz orbitals to the ground state, implying that 
orbital fluctuations play no role for the low-temperature 
physics of TiOCl. Nevertheless, the detailed comparison 
of the experimental results to our electronic structure cal- 
culations as well as to recent LDA+DMFT calculations 
show that even the normal state properties of TiOCl are 
far from being fully understood. 



II. TECHNICAL DETAILS 
A. Experimental 

Single crystals of TiOCl were prepared by chemical va- 
por transport from TiCl3 and Ti02^ The samples were 
characterized by x-ray diffraction, specific heat, magnetic 
susceptibility, and electron spin resonance measurements. 
The high quality of our crystals can be directly inferred 
from the susceptibility data displayed in Fig.^ in which a 
pronounced hysteresis at T Cl signals the first-order char- 
acter of the lower transition. 

Photoemission spectroscopy has been performed 
mainly at our homelab using He I (21.2 eV) and Al K Q 
(1486.6 eV) radiation and an Omicron EA 125 HR elec- 
tron energy analyzer. For polarization-dependent mea- 
surements we utilized an Omicron AR 65 spectrometer 
equipped with a He discharge lamp and a rotatable polar- 
izer. Additional polarization experiments with linearly 
polarized synchrotron light were performed at the SIS 
beamline of the Swiss Light Source at the Paul-Scherrer- 
Institute (Villigen, Switzerland) using a Scienta SES 100 
analyzer. The energy and angular resolution amounted 
for the homelab experiments to 60meV and ±1°, respec- 
tively, and to 80meV and ±0.2° in the SLS experiments. 
The layered structure of TiOCl facilitates easy surface 
preparation by in situ crystal cleavage. The resulting 
surfaces were clean and atomically long-range ordered 
as evidenced by x-ray induced photoemission (XPS) and 
low-energy electron diffraction (LEED), respectively. Be- 



cause TiOCl is an insulator, all photoemission data were 
taken at room temperature or above in order to min- 
imize sample charging. From systematic temperature- 
variations we found that the charging is almost negligible 
at and above « 370 K. In this situation the maximum of 
the Ti 3d-peak (see below) is located at 1.45 eV below 
the experimental chemical potential ^i exp , which corre- 
sponds to the Fermi edge position of a silver foil. Spec- 
tra measured at lower temperatures have been aligned 
accordingly. 

B. LDA+U calculations 

On the theoretical side, we have determined the elec- 
tronic structure of TiOCl by performing density func- 
tional theory (DFT) calculations in the generalized 
gradient approximation (GGA)i& and in the so-called 
LDA (GGA)+U approximation^ using the full-potential 
linearized augmented plane-wave code WIEN2kiS£ Since 
the LDA+U calculations are performed on a spin- 
polarized state, we considered a ferromagnetic alignment 
of the Ti spins (FM) as well as an antiferromagnetic ar- 
rangement of the Ti spins along the b direction (AFM). 
The latter was found to be lower in total energy compared 
to the FM stated We also performed LDA+U calcula- 
tions for the low-temperature crystal structure of TiOCl 8 
in an antiferromagnetically spin-polarized state. In all 
calculations we used RK max = 6 and 56 k irreducible 
points for Brillouin-zone integrations. The values for the 
onsite Coulomb repulsion U as well as onsite exchange Jo 
were taken to be 3.3 eV and 1 eV, respectively, 5 - which 
accounts well for the intersite chain exchange constant 
derived from the magnetic susceptibility. 

C. Dynamical density- matrix renormalisation 
group calculations 

To better account for correlation effects beyond the 
LDA+U approach we have also determined the one- 
particle spectral function of the one-dimensional (ID) 
Hubbard model at half band-filling using the dynam- 
ical density-matrix renormalisation group (DDMRG) 
methodi2ii22 The calculations were performed on 32-site 
chains with open boundaries and we kept up to 200 
density-matrix eigenstates. The broadening of the spec- 
tra is rj/t = 0.2. 

III. RESULTS AND DISCUSSION 
A. Valence density of states 

Figure [5] shows angle-integrated photoemission spec- 
tra measured with 21.2 eV (UPS) and 1486.6 eV (XPS) 
photons. Both spectra consist of two well separated re- 
gions. The small intensity differences between the UPS 
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FIG. 2: (Color online) Angle-integrated photoemission spec- 
tra and LDA+U density of states for a ferromagnetic spin- 
polarized state . The UPS and XPS spectra were measured 
at photon energies of 21.2 eV and 1486.6 eV, respectively. 



and XPS data are due to different photoemission cross 
sections. Also shown is the occupied part of the density 
of states (DOS) obtained from the LDA+U calculation 
for the room-temperature structure assuming ferromag- 
netic spin alignment (LDA+U FM ). From the comparison 
of theory and experiment we identify the low-binding en- 
ergy peak as Ti 3d-like, whereas the states between —9 eV 
and — 4eV are predominantly derived from CI 3p and O 
2p levels. The overall shape of the theoretical DOS and 
the experimental spectra are rather similar. Note, how- 
ever, that the LDA+U DOS is shifted to higher binding 
energy in order to align the theoretical and experimental 
Ti 3d peaks. With this alignment the relative separa- 
tion between Ti 3d and the ligand p states still appears 
« 0.5 eV too small. The value of 2 eV for the correlation 
gap as determined from optical spectroscopy 23 ' 24 is in 
reasonable agreement with our photoemission data and 
the LDA+U calculations (not seen here, as the figure 
displays only the occupied DOS). 

In the following we focus on the Ti 3d part near 
the chemical potential, which is shown as a blow-up 
in Fig. |3J According to the LDA+U calculations^ it 
has exclusively 2>d xy orbital character. As seen in the 
upper panel of Fig. [31 the experimental 3d spectrum 
has a similar asymmetry as the LDA+U FM DOS but 
is much broader. We emphasize that the experimental 
width is perfectly reproducible and not due to instru- 
mental broadening, which is much smaller. A comparison 
to the LDA+U DOS for antiferromagnetk spin configu- 
ration (which seems the more natural choice consider- 
ing the Bonner-Fisher-type susceptibility) yields an even 
stronger disagreement. 

We note that LDA+U accounts for the onsite Coulomb 
interaction only in a mean-field way and is thus effectively 
still a one-electron theory for static ordered systems. One 
may thus speculate that the origin for the conflicting 3d 
widths lies in pronounced electronic correlation effects 
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FIG. 3: (Color online) Ti 3d part of the angle-integrated 
photoemission spectrum (hi; = 21.2 eV) compared to the 
density of states obtained from various calculations. Up- 
per panel: LDA+U with ferromagnetic and antiferromagnetic 
spin order (U — 3.3 eV). Lower panel: ID Hubbard model, 
LDA+DMFT using iterative perturbation theory (Ref. Il7l) . 
and LDA+DMFT using quantum Monte Carlo (Ref.^. All 
spectra are normalized to same integrated area. 



and/or fluctuations of spin-Peierls or orbital nature be- 
yond the scope of LDA+U. Electronic correlation effects 
can in principle be accounted for by a combination of 
LDA and dynamical mean field theory (LDA+DMFT) or 
by suitable many-body models. The lower panel of Fig. [31 
shows a comparison of the experimental Ti 3d spectrum 
to two different LDA+DMFT calculations and the DOS 
of the half-filled ID single-band Hubbard model calcu- 
lated using the DDMRG method (using U = 3.3 eV and 
a transfer integral of t = 0.23 eV, corresponding to the 
LDA+U calculations). The LDA+DMFT calculations 
are taken from Refs. ^3 an d Hf| and used different im- 
purity solvers and basis sets. While the LDA+DMFT 
curves indeed display a broadening much closer to that 
of the experiment, none of these curves can sufficiently 
explain the shape of the experimental data. The striking 
disagreement between the two LDA+DMFT calculations 
should probably be traced back to the fact that, contrary 
to the multi-orbital QMC solver, the IPT solver is quite 
uncertain for anisotropic multi-band problems, 26 as it is 
the case for TiOCl. Also the ID Hubbard model does 
not reproduce the bandwidth and the detailed spectral 
shape of our photoemission data. 





-3 -2 -1 -3 -2 -1 




FIG. 4: (Color online) ARPES intensity plots I(k,E) and 

EDCs along the crystallographic axes a and 6, corresponding FIG - 6: (p olor online ) ARPES intensity plot 7(k, E) and 
to the XrX and YTY lines in the orthorhombic Brillouin zone. LDA+U bands along the crystallographic a-axis. 



B. Dispersions and anisotropy of the electronic 
structure 

We now turn to our angle-resolved photoemission data, 
which are shown in Fig. 0] in a broad energy range for 
the two crystallographic directions a and 6. Both data 
sets are displayed in two alternative ways: as intensity 
plots I(k, E), and as energy distribution curves (EDCs). 
As seen, there is well-pronounced dispersion particularly 
in the CI 3p/0 2p part of the spectra between — 9eV 
and — 4eV. This behavior and the clear symmetry of the 
dispersions with respect to the T point is indicative of 
the good crystal surface quality. 

Figure [3] shows an enlargement of the Ti 3d ARPES 
spectra. From a comparison of the bare EDCs one can 
easily see that the spectral changes along the &-axis are 
considerably stronger than along the a-axis. This ob- 
servation confirms the quasi-one-dimensional electronic 
(and magnetic) nature of TiOCl and corroborates the 
LDA+U prediction 5 that ID chains are formed by direct 
hopping via Ti 3d xy orbitals along the 6 direction. We 
note, however, that the dispersion along the a-axis is still 



finite as seen from a closer inspection of the fc-dependence 
of the peak maximum, which moves to slightly higher 
binding energy from T to X. The size of the a-axis disper- 
sion is directly related to the interchain coupling, which 
may play an important role for the complex spin-Peierls 
transition behavior^ A comparison of the experimental 
a-axis dispersion to that of the LDA+U FM calculation is 
displayed in Fig. Eland yields qualitative agreement con- 
cerning size and direction of the dispersion. Note again, 
however, that the LDA+U bands had to be shifted in 
energy to give a good match to the ARPES data. 

The spectral changes in the 6-axis ARPES data are 
somewhat more complex (right panel of Fig. [SJ). They 
start out at T with a single peak at w 1.5 eV below the 
chemical potential. With increasing momentum the peak 
shifts clearly towards /i exp and reaches its smallest bind- 
ing energy about halfway between T and the zone edge. 
For even larger momentum the peak rapidly drops in in- 
tensity and seems to move back to slightly higher bind- 
ing energy. At the same time a new feature appears at 
~ — 2.5eV until at the zone boundary (Y) the spectral 
shape has evolved into a broad hump. This behavior and 
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FIG. 7: (Color online) ARPES intensity plot I(k,E), LDA+U FM/AFM bands calculated for the high-temperature structure 
(U — 3.3 eV, Jo = leV), LDA+U AFM bands calculated for the dimerized low-temperature structure, and one-particle spectral 
function of the Hubbard model calculated by the DDMRG method (U = 3.3 eV, t = 0.23 eV). 



the relatively large broadening of these structures have 
been reproducibly observed on many different samples 
and is hence to be taken as intrinsic. 

Figure [7| shows the experimental 6-axis dispersion 
as ARPES intensity plot in comparison with vari- 
ous theoretical calculations. Starting with the theo- 
retical LDA+U FM dispersion calculated for the high- 
temperature (non-dimerized) structure we find pro- 
nounced disagreement with the experiment, particularly 
in the second half of the TY direction where the the- 
oretical bands continue to disperse upwards, whereas 
the experimental dispersion turns over and bends down- 
wards until at the Y-point the spectral weight distri- 
bution is strongly broadened and seems to reach even 
higher binding energy than at T. Using antiferromagnetic 
spin alignment in the LDA+U approach doubles the unit 
cell and hence results in a dispersion symmetric about 
irY (center panel of Fig. |7|). Although LDA+U AFM 
thus reproduces the dispersion maximum halfway along 
TY, it does clearly not account for the asymmetric ex- 
perimental behavior and the strong broadening towards 
the Y-point. Motivated by speculations that the unusual 
high-temperature behavior of TiOCl could be caused by 
fluctuations of the spin-Peierls order parameter, Fig. [7| 
shows LDA+U AFM dispersions for the dimerized low- 
temperature phased We note that such fluctuation ef- 
fects have been observed in charge Peierls systems well 
above the actual transition temperature. 27 However, as is 
evident from the figure, the effect of (fluctuating) dimcr- 
ization results in band doubling but is otherwise rather 
small and hence cannot explain the phenomenology of 
the ARPES data. 

The right panel of Fig. [7| finally displays the 
momentum-resolved spectral weight distribution of the 
ID single-band Hubbard model, calculated within 
DDMRG for the U and t parameters of the LDA+U cal- 
culations. We note that in this case the spectra are en- 
tirely of incoherent nature and correspond to momentum- 
dependent continua, which are structured in intensity 



due to the phase space available for decomposition of 
a real (photo-)hole into separate collective spinon and 
holon excitations. As a consequence a number of dif- 
ferent dispersive structures appear, whose detailed (non- 
quasiparticle) nature are discussed in, e.g., Refs. |3| and 
l2l Comparing the DDMRG result with our ARPES data 
we find several corresponding features such as the initial 
upward dispersion (due to spin and charge branches in 
the Hubbard modelSS), the dispersion maximum at ^TY, 
and the asymmetric shift of weight to larger binding en- 
ergy towards the Y-point. Also the overall energy width 
of the spectral weight distribution is greater than in the 
LDA+U calculations. On the other hand, the experiment 
does not show the pronounced spin-charge splitting be- 
tween T and ^Y predicted by DDMRG, nor the holon 
"shadow band"— which in the Hubbard model spectrum 
disperses downwards from T. 

In addition, we calculated the spectral functions of 
two closely related Hubbard-type models in the Mott- 
insulating phase: The extended Hubbard model which 
includes a next-neighbor Coulomb repulsion V and the 
t-t'-U model— that takes into account a next-neighbor 
hopping t' . The resulting spectral functions (not shown) 
for V/t < 2 and t' jt < 0.5, respectively, display no quali- 
tative differences compared to the simple Hubbard model 
spectral function. However, the spectral weight of the 
spin branch in the t-t'-U model is significantly reduced 
for appropriate parameters {U/t = 15.7, t'/t = 0.14) 6 . 
This may indicate why the spin-charge splitting is not 
resolved experimentally. 

Overall, the agreement is much better for our Hub- 
bard model calculations than for the LDA+U band dis- 
persions, suggesting that the experimental spectra are 
indeed dominated by electronic correlation effects. We 
note, however, that we cannot really expect the ID single- 
band Hubbard model to account fully for our spectra, as 
it completely ignores the orbital degrees of freedom. The 
crystal-field splitting between the d xy ground state and 
the excited d xz . yz states is only a few 100me\&22i and 
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hence virtual excitations into these states will be even 
more important than the double occupations on the same 
site (with energy scale U = 3.3 eV) already contained in 
the single-band model. Therefore, it would be highly 
desirable to compare our ARPES spectra to that of a 
suitable multi-orbital ID Hubbard model which however 
is still out of reach for the DDMRG. 



C. Polarization-dependent photoemission 

As already discussed in the introduction, the high- 
temperature phase of TiOCl is characterized by anoma- 
lous broadening of the phonon lines in Raman and in- 
frared spectroscopy i 11 ' 12 This indicates strong coupling 
of the lattice to the electronic (spin or orbital) degrees of 
freedom. Indeed, LDA+U calculations for frozen lattice 
distortions corresponding to the relevant phonon modes 
have shown that the orbital ground state may switch from 
3d xy to 3d xz ,y Z for sufficiently large amplitude. 6 Such a 
dynamical Jahn- Teller effect would lead to an orbitally 
mixed character of the time-averaged ground state. 

Experimental information on orbital symmetry can be 
obtained from polarization-dependent ARPES, making 
use of selection rules realized for special experimental ge- 
ometries. Assume the direction of the incident linearly 
polarized light and the emission direction of the photo- 
electrons lie within the same crystal mirror plane. Then 
with the polarization vector within and perpendicular to 
the mirror plane the ejected electrons can stem only from 
states with well-defined even or odd parity with respect 
to this plane, respectively^ If we choose the experimen- 
tal setup as sketched in Fig.[S]the dry-derived band states 
of TiOCl are even with respect to the (6,c) mirror plane 
while the d xz ^ yz states are odd. With the (6,c)-plane 
lying horizontally, photoemission from the d xy states is 
hence dipole-allowed only for horizontal light polariza- 
tion, whereas d^^-emission can be observed for vertical 
polarization only. 

In the main panel of Fig. |H1 we display the results of 
such an ARPES experiment obtained with 21.2 eV pho- 
tons at room temperature. The spectra were recorded 
at normal emission (i.e., at the T point) with horizontal 
(blue) and vertical (red) polarization, respectively. Since 
a laboratory He lamp with a rotatable polarizer was used, 
the photon flux is not changed by switching the polariza- 
tion and the spectra can be normalized to equal integra- 
tion times per channel. At first glance, one can see that 
the intensity distribution over the whole valence band 
indeed is strongly affected by polarization effects. Fo- 
cussing on the Ti 3d states reaching down to about 3.5 eV 
below /j, exp their spectral weight is significantly reduced 
but not completely suppressed upon switching the polar- 
ization from horizontal to vertical, indicating the dom- 
inance of d xy emission. The residual 1 : 4.8 weight for 
vertical polarization can be quantitatively accounted for 
by the finite degree of light polarization (« 85%), possi- 
ble small sample misalignment and the effect of thermally 
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FIG. 8: (Color online) Photoemission spectra measured at 
the F-point with horizontal and vertical light polarization, 
respectively (T = 300K). The inset shows the corresponding 
photoemission spectra measured at another photon energy 
(T = 365 K). The experimental setup is sketched in the upper 
left corner. 

activated symmetry-breaking phonons^i Similar data for 
the Ti 3c? derived part of the electronic structure using 
polarized synchrotron radiation is depicted in the inset of 
Fig- El These spectra were recorded at a photon energy of 
110 eV and a temperature T — 365 K and normalized to 
equal photon flux. Due to the higher polarization degree 
of the synchrotron radiation the residual relative weight 
for vertical polarization is even further reduced. We note 
that this contribution does not vary down to room tem- 
perature within experimental accuracy again indicating 
that phonon-induced orbital fluctuations are not impor- 
tant in this temperature range. 

Hence we have direct experimental evidence that there 
is no sizable dynamical Jahn- Teller admixture of d xz and 
d yz states to the d xy ground state at room temperature. 
This is in line with the reasoning based on recent cluster 
calculations combined with polarization-dependent opti- 
cal data—, although such an effect cannot be excluded 
for far higher temperatures. 



IV. CONCLUSIONS 

In summary, we have investigated the electronic struc- 
ture of TiOCl by means of angle-resolved photoelec- 
tron spectroscopy at and slightly above room tempera- 
ture. From those measurements we confirm the predicted 
quasi-one-dimensional nature of the low-lying energy ex- 
citations. However, the electron dispersions along the 
one-dimensional axis cannot be explained by LDA+U or 
DDMRG calculations of the one-dimensional single-band 
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Hubbard model. From the specific polarization depen- 
dence of the Ti 3d spectral weight close to the chemical 
potential sizeable phonon-induced admixtures of d xz and 
d yz derived states to the d xy ground state can be ruled 
out. Thus a scenario where the non-canonical behavior 
of the spin-Peierls transition in this compound is traced 
back to orbital fluctuations within a dynamical Jahn- 
Teller scenario does not seem plausible. 
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